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An elastic membrane that is forced to reside in a container smaller than its natural size will deform 
and, upon further volume reduction, eventually crumple. The crumpled state is characterized by 
the localization of energy in a complex network of highly deformed crescent-like regions joined by 
line ridges. In this article we study, through a combination of experiments, numerical simulations, 
and analytic approach, the emergence of localized regions of high stretching when a self-avoiding 
membrane is subject to a severe geometrical constraint. Based on our experimental observations and 
numerical results, we suggest that at moderate packing fraction, inter-layer interactions produce a 
response equivalent to that of a thicker membrane that has the shape of the deformed one. We find 
that new conical dislocations, coined satellite d-cones, appear as the deformed membrane further 
compactifies. When these satellite d-cones are born, a substantial relaxation of the mechanical 
response of the membrane is observed. Evidence is found that friction plays a key role in stabilizing 
the folded structures. 



I. INTRODUCTION 

When an elastic membrane is constrained into a small 
volume, the initially flat structure experiences Euler 
buckling [TJ [2] . As the accessible volume decreases, the 
membrane goes through a second transition generating 
a region where the strain localizes and a crescent sin- 
gularity is born. This new configuration minimizes the 
total energy at the expense of the emergence of this pivot 
point. Beyond the crescent, the membrane forms a struc- 
ture known as developable-cone (d-cone) [3-6 . Subse- 
quently, as the available space is further reduced, com- 
plex patterns appear, ridges that join these localized ob- 
jects forming a network in which energy localizes [7HT2]. 

Besides bending and stretching, connectivity and self- 
avoidance are also key ingredients in crumpling [TTJ [13] , 
both determining the evolution and the final fate of the 
folding process. On the one hand, the connectivity of a 
sheet forces the emergence of the above-mentioned net- 
work of highly stretched regions when volume restric- 
tions are imposed. On the other, self-avoidance plays a 
key role in the mechanical response of the crushed mem- 
brane since perfectly elastic membranes that incorporate 
excluded volume energy are much stiffer than phantom 
ones [11 . This suggests that it may be possible to ratio- 
nalize a self- avoiding membrane as one effectively thicker 
than its phantom version. 

In this paper, we study a model system that allows us 
to make a natural link between the folding mechanism 
and the birth of new crescent singularities. The exper- 
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imental setup consists of a circular membrane of radius 
i?, that under the action of a point-like traction load 
F directed along the vertical direction, is forced to pass 
through a hole of radius R^. The plane of the hole divides 
the space into two domains z > and z < 0, as shown 
in FigJI] We find that the mechanical response (force 
versus displacement plot, F-z, see Figj3| exhibits a se- 
ries of hardening-softening events as well as catastrophic 
events at large confinement. A direct relation between 
the emergence of these prominent non-linearities, the fric- 
tional collision of folds, and the birth of new singularities 
is reported and discussed. 

The paper is organized as follows. Section II is a de- 
scription of the experimental methods and results are 
reported in Section III. Section IV is devoted to numer- 
ical simulations as well as to the analytical results valid 
in certain regimes of the force versus displacement curve. 
Finally, Section V contains concluding remarks. 



II. MEASUREMENT OF THE MECHANICAL 
RESPONSE 

The membranes used in all our experiments are made 
of acetate. Their thickness is t ~ 0.10 mm and their 
radius R = 100 mm. Their bending stiffness was ob- 
tained by measuring the slope at the early stage of the 
experimental mechanical response curve of the membrane 
shown in Fig (3) B = (4.2 ± 0.1) x 10~ 4 N m (details of 
the scaling in this stage are presented in Section IV). 
We performed a tensile test (using the Instron 5544 
system) to determine the Young modulus of the mem- 
branes, Y = (3.30 ± 0.02) x 10 9 N/ra 2 . The Poisson 
ratio v ~ 0.35, was then computed using the formula 
Ft 3 /12(l - v 2 ) = B. The radius of the hole is fixed at 
Rh = 23.5 mm. 

To pull a membrane by its center, we punched a small 
hole at the center with a surgical needle and put a thin, 
stiff wire through the hole. At one end of the wire we at- 




Figure 1: (Color online) Typical conformations that appear 
when pulling an almost inextensible circular membrane by 
its center through a hole, (a), Low load regime in which a 
d-cone shape is observed. This is an intermediate state of 
the first regime, (b), Opposite sides of the initially flat circle 
touch and self avoidance becomes crucial. This event marks 
the onset of the second regime, (c), and (d), are lateral and 
top views of high- load regimes. The main signature of the 
second regime is shown in (c) where the mirror symmetry is 
broken. The satellite d-cone is seen in (d). The onset of the 
satellite d-cone breaks the translational symmetry along the 
conical generators. This event marks the onset of the third 
regime. In (a) the relevant variables to describe this system 
are shown: R is the radius of the membrane, Rh is the radius 
of the hole and F represents the load. 



tached a copper ball of diameter ~ 2 mm. Then, the wire 
was pulled in order to produce a contact region between 
the copper ball and the membrane. Afterward we painted 
that region with a small amount of silicone glue to restore 
the membrane connectivity. We connected the free end 
of the wire to a commercial digital dynamometer (Ultra- 
Sport) and a home-made system of screws that allowed a 
controlled displacement in the vertical direction and the 
simultaneous measurement of the reaction force exerted 
by the membrane. Special care was taken to make the 
size of the perturbed region much smaller than the typi- 
cal size of the core of the elastic singularity, which for the 
present case was r c ~ 5 mm. We pulled the membrane 
down in the vertical direction i, at a rate of 0.015 mm/s, 
a factor 30 slower than in a previous work [14] . This is an 
important point since the mechanical response depends 
on how fast the membrane is pulled because of friction. 
This dependence becomes stronger when folds collide. In 
all experiments the vertical displacement of the tip of the 
membrane was measured directly with a digital microme- 
ter Mitutoyo. We observed that each time the membrane 
was pulled a distance Az, it took at least 2 seconds for 
the sliding process to stop. For this reason, every time 
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Figure 2: (Color online) Typical conformations that appear 
when pulling a conical surface by its center with a force F. 
(a), Low load regime in which the membrane is strained, (b), 
At a critical force a satellite d-cone appears, (c), For larger 
displacement in the z direction a hinged structure appears, 
(d), One side of the cone touches the opposite one, resembling 
the geometry of Fig[TJb). 



we pulled the membrane, we waited approximately 20 
seconds until the sheets had completely relaxed. This 
interval was long enough for the reaction force measured 
by the dynamometer to reach a constant value, confirm- 
ing that the sheet had reached a stationary state. As 
the membrane progressed along the vertical direction, it 
was more difficult to pull it and the force needed to drag 
it down became larger and larger. But this process was 
non-monotonic: every time two folds collide and slide 
over each other, or a catastrophic event such as when a 
satellite d-cone is born, the force needed to pull and hold 
the membrane drops dramatically, indicating an instabil- 
ity of the packing structure of the membrane. 

To obtain data with reliable statistics, we repeated 
each experiment 10 times (and each time we used a 
new membrane of the same physical characteristics as 
described above). We studied regular, lubricated, and 
sticky membranes in order to examine the role played by 
friction. To obtain membranes with different friction co- 
efficients but the same mechanical properties, we coated 
the regular acetate membranes. The reduced friction case 
was realized by applying a thin layer of graphite powder. 
We applied the graphite with a clean brush and remove 
the loose powder. In the opposite limit, to increase the 
friction we applied a thin film of 3M spray glue, 5 min- 
utes before the experiment was started. For obtaining as 
uniform as possible films, we sprayed the glue to the sam- 
ples from 1 meter apart during 5 seconds. The change 
in the thickness of the samples due to the glue coating 
was measured using a micrometer, being in average of 
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Figure 3: (Color online) Experimental mechanical response 
performed in initially flat membranes: Force vs. displacement 
curves of lubricated (green open circles), normal (red trian- 
gles pointing up) and sticky (blue triangles pointing down) 
membranes for increasing confinement. The first regime is 
realized from z — to z ~ 25 mm, in the three types of mem- 
branes. The leftmost (blue) arrow located at z ~ 33 mm, 
shows the first time that two folds collide in the three types 
of membranes. The major catastrophic event corresponds to 
the development of a satellite d-cone (denoted by the middle 
and rightmost (black) arrows which also denotes the onset of 
the third regime), which occurs earliest for sticky membranes 
(zd-cone ~ 40 mm) and latest for lubricated (zd-cone ~ 60 
mm for regular membranes and the satellite d-cone for lu- 
bricated membranes is not shown in this plot since it occurs 
for larger z). The left upper inset shows a schematic rep- 
resentation of the collision of two folds due to the pushing 
force F/2 at each edge. A/2 denotes the displacement of 
each of the edges of the rod from its initial position along 
the boundary. The lower inset shows a zoom-in look at the 
low load regime [6] along with the best fit obtained using 
F — zB/Rh \og(R/r c ) where B is the free parameter reported 
in the text. 



9 percent. The change in thickness due to the graphite 
powder was negligible. The results of the tensile tests 
(performed with the Instron 5544 system) applied to the 
coated membranes yielded a change of the Young's mod- 
ulus due to the coating of 0.3 percent. The change in the 
bending stiffness was negligible as well, as apparent from 
FigJ3] The mechanical response for these membranes are 
shown in Figj3]with error bars deduced from variations 
of the experimental data on repeated runs. 

To examine the deformations that appear when an 
already compacted structure with conical symmetry is 
pulled, we performed a similar experiment with mem- 
branes whose initial geometry is conical, Figj2j Each cone 
was cut out of a circular membrane with a deficit angle 
of 135°. We glued opposite sides of the cut with epoxy 
and used a hair dryer to release internal tensions in the 
material. We repeated each experiment 10 times as well 
for the same deficit angle. We did observe that a d-cone 
never appeared on the suture line most because along the 
suture the bending stiffness is larger than in the rest of 



the structure. The mechanical response is shown in Fig(4] 
A more detailed description of the geometry and evolu- 
tion of these deformations will be presented elsewhere 



III. EXPERIMENTAL RESULTS 
A. Compaction of the membrane 

Experimentally we found three qualitatively different 
regimes of the mechanical response of the membrane. 
During the initial stage the membrane deforms elasti- 
cally until the pulling force F reaches a critical value F c 
at which the first d-cone singularity appears [4j [6] . As we 
keep pulling, the single-valued description used in previ- 
ous work [4] breaks down. The d-cone structure evolves 
until it touches itself, producing a contact interaction be- 
tween the two opposite radial lines of the membrane due 
to self avoidance. This new structure is stiffer and, if ob- 
served from the top (See Fig{l]3), contains three regions 
separated by the contact lines. The first contact between 
the two folds marks the end of the first regime and the 
birth of the second and occurs for z ~ 25 mm (shown in 
Fig J3] and Fig(l]3). Once in the second regime, the fric- 
tional interaction of two folds at z ~ 33 mm, produces 
a large mechanical response as denoted in Fig (3] by the 
leftmost (blue) arrow. 

As the vertical displacement is increased at a constant 
rate beyond z ~ 33 mm, the main central fold shown in 
Fig{l]3 is pushed against the opposite side of the mem- 
brane and eventually slops when friction of the two exter- 
nal folds can not stabilize the system any longer. After 
that, a main spiral develops and breaks the mirror sym- 
metry of a simple d-cone. The break of mirror symmetry 
is the signature of the second regime (see FigfTJ^). How- 
ever, in this process the translational symmetry along 
the conical generator is still preserved, particularly out- 
side the core region near the membrane tip. As we keep 
pulling along z, more folds are produced and forced to 
touch each other, the whole sheet becomes highly rigid, 
and due to self avoidance, it eventually shows an en- 
hanced mechanical response, a peak signaling the fact 
that the membrane is stuck. At this point, if we increase 
the load further, a catastrophic event will occur: the for- 
mation of what we dubbed a satellite d-cone that breaks 
the local translational symmetry along the conical gen- 
erators, as shown in FigjTJd) and Fig{2|d). This event is 
characterized by a sudden drop in the force-displacement 
curve, indicated by the middle and rightmost (black) ar- 
rows in Fig 3] (zd- scone depends on friction, ^d-scone ~ 40 
mm for sticky membranes ^d-scone ~ 60 for regular mem- 
branes) which marks the end of the second regime and 
the onset of the third one. The rearrangement of folds 
allows the system to keep evolving into more compact 
states as the membrane is being pulled along z. We em- 
phasize that the breaking of the translational symmetry 
along the conical generators makes this effect intrinsically 
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Figure 4: Mechanical response for initially conical surfaces: 
Force vs. displacement curve. The conical generator length 
is R — 100 mm and the deficit angle is — 135°. The actual 
membrane and deformation is shown in Fig[2] The arrow 
indicates the onset of a satellite d-cone. The inset shows 
a zoom-in look to the mechanical response before the first 
satellite d-cone appears (open circles) along with the best fit 
(line). In this region F ~ 2 ,i-50±o.02 > Error bars are smaller 
than the size of the data points and thus not shown here. 

three dimensional, and any attempt to model this soft- 
ening with a two dimensional model is doomed to fail. 
Indeed, a study based on a two-dimensional description 
[T4] suggested that at large compactations divergent me- 
chanical responses are expected. 

The emergence of a satellite d-cone is the striking sig- 
nature of the third regime. This is the simplest real- 
ization of an energy cascade proposed as the soul of 
crumpling a long time ago and a beautiful example of 
a nested generation of energy condensation [6j [12] . The 
three curves shown in Fig(3] clearly show that friction 
does play an important role in the mechanical response 
of our system. 

Friction in two dimensional packing of wires was shown 
to have a strong impact in the packing morphology [16], 
showing the so-called classical, spiral, and plastic mor- 
phologies. However, in two dimensions crumpled wires 
do not have the freedom to release energy by deform- 
ing into another dimension as in the case discussed here. 
In two dimensional systems, a divergent mechanical re- 
sponse close to a non-universal (material-dependent) crit- 
ical packing ratio is generically observed, which for cer- 
tain parameters was found to be <\> c ~ 0.46 [16]. Our case 
is different and does have a third dimension that allows 
new deformation modes. One example is the satellite d- 
cone observed here. We further found that the vertical 
displacement at which the satellite d-cone forms depends 
on friction. The middle and rightmost (black) arrows 
in Figj3] indicate that the satellite d-cone appears earlier 
when the friction is larger. This can be qualitatively un- 
derstood in the following way: larger friction means it is 
harder for touching layers to slide over each other, which 
makes the membrane effectively more rigid. Therefore, 



for membranes with a larger friction coefficient, the com- 
paction preserving the translational symmetry along the 
conical direction stops at smaller vertical displacements. 
As the displacement is further increased, the satellite d- 
cone appears to allow the membrane to further compact. 

All our experimental and simulation results indicate 
that the satellite d-cone is the basic deformation mech- 
anism of a conical surface, which is effectively the ge- 
ometry of the flat membrane when it is in a compacted 
state with translational symmetry. To provide more evi- 
dence of this statement, we performed an experiment of 
a membrane whose initial geometry was conical. This ex- 
periment (see Figsj2|4| shows that for large confinement, 
once the folded membrane is stuck at a certain friction 
dependent position, the system starts to stretch until it 
reaches a critical force F c at which energy focuses in a lo- 
calized region. This region is at a certain distance away 
from the membrane tip and is the core of the satellite 
d-cone. 

The mechanical response of the cone is shown in Fig(4] 
The first stage is characterized by a dragging force that 
grows as a power law F ~ z q , where q = 1.50 ± 0.02, as 
obtained from the best fit of the F — z curve shown in 
the inset of Fig(4] In this regime the overall symmetry 
of the cone remains unchanged, therefore the measured 
force can only be due to the stretching of the membrane. 
After a critical load F c has been reached, or equivalently, 
after a critical amount of energy has been stored in the 
stretching mode, the conical shape becomes unstable and 
a satellite d-cone is born, as shown in Figj2}} and as in- 
dicated by the arrow in Figj4] The energy transfer be- 
tween stretching and bending modes is reflected in the 
F-z curve as a sudden drop in F, similar to what is ob- 
served in Fig |3] for a flat membrane in a compacted state. 
After the satellite d-cone is formed, the system evolves in 
a similar way to the evolution of the initially flat mem- 
brane at intermediate packing fraction. In this regime 
the F-z curve goes down until opposite sides of the cone 
touch each other. 



B. The role of friction 

Figj3] shows that at small vertical displacements, fric- 
tion and self-avoidance play no role because different 
parts of the membrane do not touch yet. In this region, 
all membranes with different frictional treatments exhibit 
the same mechanical response and the force-displacement 
curves collapse. However, we do not anticipate any pos- 
sible collapse of the F-z curves beyond the first fold col- 
lision. In this region, the different friction properties of 
the membranes affect the collision and sliding of folds, 
leading to friction-dependent mechanical responses, like 
the non-universal behavior found in the compactation of 
wires in 2D [16]. 

Generally, for larger friction wider folds are observed 
before they slip. This observation allowed us to further 
unveil the role of friction quantitatively. We estimated its 
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contribution by computing the friction coefficient in two 
ways. In the first approach we analyzed a fold collision at 
low packing fraction (which is representative of the pro- 
cesses happening at more compact configurations). An 
example is shown in Fig [I}}, where two folds touch each 
other. This state will last as long as the forces in the di- 
rection perpendicular to the contact line do not overcome 
the frictional forces. As the pushing force increases due 
to increasing confinement, both folds become wider and 
eventually buckle as if they were a single rod, forming 
an effective Euler elastica (see schematic representation 
shown in the inset of Figj3J) . At this stage the contact 
point will work as a joint and the folds will move with 
respect to each other: the reaction force will increase un- 
til they slip. We make use of this description to estimate 
the stability of the collision of two folds in the experi- 
ment. The only force stabilizing the system is friction, 
/iN, in which \i is the static coefficient of friction and N 
is the normal force. By virtually cutting an Euler elastica 
into two halves, we compute the bending force needed to 
keep the system together with the virtual work technique. 
Equating these two forces yields 



2(1 -v 2 ) 



(1) 



where L is the size of the effective rod, v ~ 1/3, its Pois- 
son ratio, and A is the displacement of the edges of the 
rod from its initial position along the boundary. Experi- 
mentally we found A/L ~ 1/15. Thus \i can be estimated 
using Eq{T] In the second approach, we obtained \i using 
a traditional method: An incline and a block that is free 
to slide as soon as the incline reaches a critical angle C . 
The friction coefficient in this case is well known to be 
given by 



tan (0 C 



(2) 



The coefficients of friction obtained by these two inde- 
pendent methods agree well within experimental error 
/i = 0.45 ± 0.05. This indicates that the set of most 
prominent bumps observed in Figj3j at low packing frac- 
tion are a consequence of the frictional collision between 
folds and that the height and position of these events in 
the F-z curve becomes material dependent since they are 
determined by friction. 

In our experiments using membranes lubricated with 
graphite powder, the friction is very low and the onset 
of the satellite d-cone occurs for a larger load than in 
the case of the regular membranes. In other words, the 
satellite d-cone occurs later in the F-z curve. To study 
the opposite limit where friction is extremely large, we 
deposited a thin film of a 3M glue on top of the acetate 
membrane. This treatment makes the folds stick once 
they touch each other and they remain stuck unless the 
adhesive layer fails. But the failure never occurred in 
our experiments. In this extreme case, the membrane 
can only invaginate from the configuration in which the 
two main folds shown in Figjl}} are glued together. We 



discovered that the effect of this strong attractive con- 
tact interaction is to produce a cone like the one shown 
in FigfTJ) earlier than for the regular membranes. Conse- 
quently, the onset of the satellite d-cone occurs earlier in 
the F-z curve (see middle (black) arrow in Figj3|. 



IV. MOLECULAR DYNAMICS SIMULATIONS 

To understand the compaction of the membrane in 
more detail, we modeled a membrane dragged through a 
small hole using molecular dynamics (MD) simulations. 
This allows us to explore and analyze not only the me- 
chanical response of the membrane, but also its energetics 
and compacted geometry. 

We used the computational model of an elastic sheet 
constructed by Seung and Nelson [17 . The sheet was 
modeled as a triangular lattice with neighboring vertices 
bonded by springs. The harmonic potential of each bond 
is 



U hoil d(r) = ^k(r - d) 2 , 



(3) 



where k is the spring constant and d the equilibrium bond 
length. We express all simulation parameters and results 
in terms of the Lennard- Jones (LJ) units based on the 
LJ potential defined as 



VLj(r) =4e LJ 



V r 



12 



(4) 



where 6lj and <tlj are the characteristic bending energy 
and length. Their units are denoted as e and a, respec- 
tively. LJ units are based on e, <r, and the mass m of the 
vertices. For example, the unit of time is r = ^mcr 2 / e. 
For the bond potential in Eq.[3]we used k = 5.7 x 10 5 e/a 2 
and d = 2 1//6 a. The bending energy of the sheet is given 

by 



E/bend(<W = •/(! + < 



(5) 



where J is the characteristic energy of bending and <p is 
the dihedral angle between two adjacent triangles that 
share a bond. For simulations reported here, J = 2.85 x 
10 3 e. From k and J, one can derive the stretching and 
bending moduli of the sheet by identifying the low energy 
limit of this lattice model [1 7 . In our case the bending 
modulus is B = 2.46 x 10 3 e and the 2-D Young's modulus 
Y — 3.29 x 10 6 e/a 3 . The equivalent sheet thickness is 
given by t = ^SJ/k = 0.2cr. 

The self- avoidance of the sheet is imposed by adding an 
extra interaction between all non-bonded vertices. This 
interaction is described by an LJ potential with 6lj = 
lOOe and <tlj = cr. The potential is truncated at r = 
2 1 / 6 <t to make it purely repulsive. The same potential 
is used to impose the interaction between the sheet and 
the confining hole. The later is implemented as a ring of 
vertices uniformly distributed on a circle with radius 
and spacing 0.05d 




Figure 5: (Color online) Configurations from molecular dynamics simulation, for Rh = 25cr, R = 290. 7cr, k = 5.7 x 10 5 e/cr 2 
and J — 2.85 x 10 3 e. The (a.i) plots (i — 1,6) show the 3-D shape at different packing stages: (1), The appearance of the 
first d-Cone. (2), The first touch of the two branches of the d-Cone. (3), The two branches start to slide over each other. (4), 
A highly packed state with many folds. (5), The appearance of the satellite d-cone. (6), Before the membrane escapes the 
confining hole. The (b.i) plots show the cross-section of the membrane in the plane of the confining hole. Figures (c,i), (d.i), 
and (e.i), Are the density plots of the stretching, bending, and self-avoidance (excluded volume) energy, respectively, projected 
back to the plane of the undeformed membrane. 



The MD simulations were performed using the Large- 
scale Atomic/Molecular Massively Parallel Simulator 
(LAMMPS) developed at Sandia National Laboratories. 
A velo city- Ver let algorithm with a time step dt = 5 x 
10 _4 r was used to integrate the equations of motion. We 
performed two types of simulations. In one, a gradually 
increasing force was applied to the tip (the central vertex) 
of the sheet to pull it through the hole. At a given force, 
MD was run for a certain period (usually lOOOr ~ 5000r) 
and then an energy minimization based on the Polak- 
Ribiere version of the conjugate gradient algorithm was 
performed. In this type of simulation, the energy of the 
sheet is minimized at each given force. In other simula- 
tions, the sheet was dragged through the hole by pulling 
its tip at a constant velocity v z along the vertical direc- 
tion. The resistance force on the tip from the other part 
of the sheet was measured. We went to very low v z to 



reduce the inertial effect. In this case the measured force 
was found to be equal to the force exerted on the sheet by 
the confining hole within statistical errors. This indicate 
that the inertial effect is negligible. 

The sheet is kept at a very low temperature T = 
10 -4 e/&B through a Langevin thermostat. The Langevin 
damping rate T is O.lr -1 for low v z . For high v z , higher 
damping rates have to be used. Tests with different T 
showed that results reported here are not sensitive to the 
damping as long as the inertial effect is negligible. 



A. Numerical Results 

In Fig(5]we summarize the evolution of the geometry 
of an in-silico membrane that is being pulled through a 
hole. The results shown were obtained by pulling the 
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Figure 6: (Color online) Mechanical response of an in- 
silico membrane from MD simulations: Force vs. deforma- 
tion curves at three different pulling velocities, v = 0.002 
(continuos red line), v = 0.005 (yellow dashed line) and 
v = 0.01 (black dotted line), for R h = 25a, R = 290.7a, 
k = 5.7 x 10 5 e/a 2 and J = 2.85 x 10 3 e. The solid (black) 
line is from a simulation where the dragging force is increased 
with small steps and the displacement of the membrane tip 
is measured. In this latter simulation MD is followed by en- 
ergy minimization to allow the membrane to fully relax at any 
given dragging force. 



central vertex of the sheet at a constant speed. The con- 
figurations of the membrane shown in row a indicate that 
the model does capture all the regimes found in the ex- 
periments: The formation of the d-cone, collision of the 
main central fold, FigjE] (a/2), and the subsequent for- 
mation of a spiral-like structure with folds touching and 
sliding. Finally the formation of the satellite d-cone is 
clear from FigsjE] (a — b/5 — 6). 

To qualitatively analyze how energy is transferred 
between different modes (bending, stretching, self- 
avoidance) we also show density plots of the different 
energy densities. It is clear from Fig (5] (c/4 and d/4) 
and Fig{5](c/5 and d/5), that before the satellite d-cone 
is formed stretching concentrates around the tip and at 
the confining ring. After the satellite d-cone is formed, 
stretching goes down and the membrane folds. Energy 
is transferred from the stretching mode to the bending 
mode (see bright scars in Fig{5] (d/5). The results in 
Fig. [5] indicate that the distributions of the bending and 
stretching energies are very heterogeneous and they are 
strongly localized near the tip of the membrane, the apex 
of the satellite d-cone, and the ridges along which the 
membrane is bent when the packing ratio is very high. 

Overall, the mechanical responses of the in-silico mem- 
brane shown in Figj6] are qualitatively similar to those 
observed in our experiments. The membrane initially 
deforms through stretching but the conical symmetry is 
preserved. The conical symmetry breaks down when the 
d-cone appears, though there is still a translational sym- 
metry along the conical generators, at least outside the 
core of the d-cone. A slight drop in the dragging force ac- 



companies the formation of the d-cone. Then the front of 
the d-cone moves inward and its two branches approach 
and collide. After this collision, the front of the d-cone 
keeps moving towards and finally collides with the oppo- 
site side of the membrane. After this point the two con- 
tacting branches of the d-cone start to flatten in response 
to the increasing confinement. At certain dragging force 
the two branches start to slide, inducing a clear drop in 
the measured dragging force on the membrane. This cor- 
responds to the first major bump in the F-z curve. When 
the confinement is further increased, more branches and 
folds are born, collide, and slide over each other. Each 
sliding event accompanies a drop in F, appearing as a 
bump in the F-z curve. But overall F increases as the 
deformation grows. It should be emphasized that the 
translational symmetry along the conical generators is 
still preserved in this stage. 

When the displacement of the membrane tip exceeds a 
certain amount, or equivalently when the confinement is 
strong enough, the translational symmetry breaks down 
with the appearance of a satellite d-cone. The apex of 
this second d-cone is away from the membrane tip. It 
appears in the form where many folds buckle and bend 
towards the interior of the compacted membrane. At this 
catastrophic event, the dragging force F drops dramat- 
ically by an order of magnitude, shown as a steep drop 
in the F-z curve. After this event, F starts to increase 
again as the confinement is further increased. 

Though quantitative agreement is not expected, there 
is also an obvious qualitative difference in the mechan- 
ical responses of the membrane between numerical and 
experimental results. For example, in the second regime 
after the birth of the d-cone and the appearance of the 
satellite d-cone, more bumps in F are observed in the 
numerical results (three bumps marked with three blue 
arrows in the simulations of FigjT^i, instead of the one 
observed in experiments). These bumps correspond to 
the consecutive collision and sliding of folds. The reason 
for this difference is that in-silico membranes are softer 
and less frictional. Figj3] shows that the less frictional 
membrane, the one lubricated with graphite, shows more 
bumps than the others, closely resembling the numerical 
results. However, it should be emphasized that due to 
the roughness of the iso-potential lines that describe the 
surface of the in-silico membrane, there is an effective 
friction, even though we did not put any explicit friction 
law in the simulation. The effective friction is an analog 
of the tangential force that would appear when we try to 
slide two egg boxes that are on top of each other. Each 
atomic stick-slip process appears as a sudden bump in 
the F-z curve. This is also the reason that the numer- 
ical results are noisy, albeit quite reproducible for each 
dragging velocity. 

Several experiments on crumpling were performed by 
compressing sheets with increasingly heavier weights [18] . 
In our set-up, this corresponds to applying a dragging 
force on the membrane that increases gradually. Results 
showed that this approach gives results different from 
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those from experiments where the confinement or defor- 
mation of the sheet is controlled. It should be empha- 
sized that experiments performed with increasing forces 
do not provide detailed information about the relaxation 
of the membrane, while such information is available from 
experiments performed with deformation controlled and 
mechanical response measured. 

Simulations also allow us to vary the dragging force 
F and measure the vertical displacement z of the mem- 
brane tip. The latter is one signature of the mechanical 
response of the membrane. The solid (black) curve in 
Figj6] shows one F-z curve generated in such a simula- 
tion. In the simulation the dragging force is increased in a 
stepwise manner. In most cases, the membrane reached a 
steady crumpled state quickly and z was easily obtained. 
However, at certain values of F, the membranes took a 
long time to relax. Typically, each relaxation corresponds 
to the transition from one locally stable compact config- 
uration to another. For example, when the load reached 
the critical value at which the d-Cone appeared, and was 
kept at this point, the fold grew and moved towards the 
opposite side of the membrane, and its two branches ap- 
proached each other. This process ends when the two 
branches touched and the front of the fold contacted the 
opposite side of the membrane. The membrane entered 
a new configuration, stabilized by the effective friction. 
This transition corresponds to the first plateau in the 
F-z curve as shown in Figj6] Successive plateaus were 
observed at the occurrence of strong reorganization of 
folds, including the birth of new folds or abrupt interfa- 
cial sliding of contacting folds. The results here make it 
clear that the full F-z curve can only be obtained when 
z, not F, is controlled. The underlying reason is simply 
that z is not a single-valued function of F due to the 
hardening-softening of the membrane. 

B. D-cone formation 

In real experiments it is extremely difficult to measure 
the initial deformation of the membrane. Intuition tells 
us that for any real membrane at low loads it should 
stretch, and then, once the stored energy is sufficient to 
produce a global change, a d-cone will be formed. This 
regime is readily explored for in-silico membranes. Re- 
sults show that the F-z curve grows in a quasi-parabolic 
fashion, F ~ z 1,33±0,02 . This indicates that pure bending 
is not the relevant mode at extremely low loads. Oth- 
erwise, F ~ z as predicted for a perfectly inextensible 
membrane [6]. 

In Figj7]3, the F-z curve at the early stage of the fold- 
ing process is shown. Here numerical data for four dif- 
ferent Rh are plotted. Results show that the force has a 
universal functional form given by 

F = f(z)^r (6) 
This is a clear signature that at extremely low loads the 
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Figure 7: (Color online) Mechanical response of the mem- 
brane for increasing confinement: (a), Points (green) corre- 
spond to one set of numerical data (N), showing many stick- 
slip events (from the left first, second and third arrows) as well 
as a huge relaxation for z ~ 175d denoted by the rightmost 
(black) arrow. The darker solid (black) line is from the aver- 
aged data (Av) and the monotone (red) line (A) is the best 
fit obtained from the model described in the text. The inset 
shows a schematic representation of the geometry used for the 
analytic calculation of the mechanical response at moderate 
packing, where z is the vertical distance from the tip, Rh is 
the radius of the hole and r is the polar radius of the cone at 
a vertical distance x from the tip. (b), In the low load regime 
the F — z curves collapse into a single universal function f(z) 
defined by Eq.(|6]). Rh denotes radius of the hole. 



membrane first stretches, and the d-cone only appears 
when a critical load F c is reached at z ~ 5d. After the 
occurrence of the d-cone, the main energy contribution 
comes from the bending of the membrane, and the scal- 
ing shown in Eq(6] breaks down. This is seen from the 
discontinuity of the mechanical response at z ~ bd that 
marks the onset of the first d-cone. The correct scaling 
for the force just after the first d-cone is formed is well 
known, F ~ zB/R\ log(i?/r c ), where B is the bending 
stiffness and r c = 5 mm is the core radius which is of the 
order of the crescent singularity. This scaling works for 
small deformations (small confinement) z/Rh < 1/4 [6]. 
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C. Mechanical response at moderate packing. 



The bending energy from this spiral region is: 



We used a phenomenological route to explain F-z at 
moderate confinement after the first d-cone is formed and 
before the emergence of the satellite d-cone ( 5d < z < 
175d in Figj7^i) in terms of the bending energy of the 
system only. A similar approach was used by Boue et al. 

We consider the bending energy in this regime as the 
sum of the bending energy of two regions that fold, qual- 
itatively, in a different way (look at the inset of Fig{7^i 
for clarity). The first contribution to the bending energy 
is from the isolated big fold living in the unoccupied core 
region. The second region contributing to the bending 
energy is the sector of the membrane far from the center 
of the cone. 

The length of the enclosed inextensible membrane at 
the hole is: 



L determines the maximum polar angle: 



tax ~ 27TW —j + 1 



(7) 



(8) 



Since the bending energy of the core region scales with 
the square of its curvature, we have that: 



EC ~ / / RoLdsdz 



(9) 



The characteristic radius for the fold living in the core 
region can be estimated as: 



R c 



X 



(10) 



where r is the polar radius of the cone at a vertical dis- 
tance x from the tip and t is the thickness of the mem- 
brane. Placing this back into Eqj9]we get the final form 
for the bending energy of the core: 



7 b c °- d ~ irBty/QTFW+F) x 
(R ~ r c ) 



log 1 



(ii) 



where R is the radius of the membrane and £ = zj . In 
the limit when t —> 0, the logarithmic contribution ap- 
proaches log (R/r c ). This simple model yields the signa- 
ture of a divergent response when r c = t(l+£ 2 ), evidenc- 
ing that there is a critical distance z at which becomes ex- 
tremely expensive to deform the core region. The critical 
pulling distance is given by: z c = Rh\/r c /t — 1 ~ 165d 
which compares well with the numerical data (z = 175d 
see Figj7^L, black arrow). 

The second region is illustrated by the cross-section cut 
of the membrane shown in Fig(2] We model this region 
as an archimedian spiral. Its radius is 



^spiral ~ \]r 2 -2t^T^ + 



to 

2^ 



(12) 



log 



(l + r c /*) + £ 2 



(13) 



Finally, the total energy (approximated as the bending 
energy only) at moderate packing, U^^ tlc ~ c^bend + 
^bencf wnere a is a numerical constant to be determined. 
Estimating r c = Id from the geometry of the deformed 
sheet from the numerical simulation (the size of the cres- 
cent in Figj5| and considering the constant a and the 
bending stiffness B as a free parameters in EqsjTT] and 
[l3j we have performed a numerical best fit of our ana- 
lytic expression U^^[ tlc and compared it with the nu- 
merical bending energy shown in the top curve of FigjH^i. 
From this operation we obtained B = 2.463 x 10 3 e and 
a = 3.26. Using these values in U^^[ tlc , we computed 

the total force in this regime, ^ t o^!|j^ tlc , which is shown 
in the monotone (red) curve in Figwa. Even though we 
have completely neglected self-avoidance and stretching, 
we obtained qualitative agreement in the mechanical re- 
sponse between our analytic approach and the numeric 
results, up to the point at which a satellite d-cone devel- 
ops. This agreement indicates that in this regime the de- 
formation of the membrane is dominated by the bending 
mode. This observation is further strengthened by de- 
composing the total energy into bending, stretching, and 
self- avoidance contributions, Figj8^i, which shows that 
bending energy is larger than the other two by an or- 
der of magnitude. 



D. First and second generation d-cones 

The catastrophic relaxation event in the F-z curve 
shown by the rightmost (black) arrow in Figj7^i is the sig- 
nature of the birth of a second generation d-cone in the 
real and in-silico membranes. The mechanical response 
of the system after this new singularity was born can be 
rationalized in terms of a renormalized bending stiffness 
of the sheet which now resembles a cone that is thicker 
than the original membrane. Indeed, since bending en- 
ergy of a d-cone scales as ZTbend ~ Bz 2 /R% log (R/r c ) for 
small deformations, we expect that the ratio between the 
response of the first d-cone and the second d-cone will be 
proportional to (t e ff/t) 3 , where t e ff = nt and n is the 
number of folded layers, apart from geometric factors due 
to the change in the contact angle. From numerical sim- 
ulations results shown in FigjFjwe see that n ~ 4 at the 
onset of the second d-cone. By comparing, on the other 
hand, the slope of the numeric Ub-z curve (Figj8^i top 
curve) we get n ~ (J^-) 1 ^ 3 ~ 3.2. The agreement is rea- 
sonable, considering that the four layers are not glued 
together. Note that membranes with very small coeffi- 
cient of friction will deviate from this simple estimate, 
since in that case the folded layers are not going to show 
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Figure 8: (Color online) Development of different inter- 
actions for increasing confinement. (a), Stretching (mid- 
dle curve, red triangles pointing up), bending (top curve, 
green circles), and self-avoidance (bottom curve, blue trian- 
gles pointing down) energy as a function of vertical displace- 
ment, (b), Individual contributions to the total force arising 
from bending (top, green circles), stretching (middle, red tri- 
angles pointing up), and self- avoidance (bottom, blue trian- 
gles pointing down) interactions. Transitions between differ- 
ent regimes at z ~ 5d and z ~ 185d indicate the formation of 
the first d-cone and the satellite d-cone. 

a cooperative behavior but most likely will slide respect 
to each other. 

Using our argument of the renormalized stiffness con- 
stant, we turn to compute the number of folded layers for 
which the satellite d-cone would be created. Examining 
FigjHfD we see that bending and stretching forces become 
comparable at the onset of the new d-cone. The reason is 
that once the membrane has become thick enough and is 
stuck due to the geometrical constraint, it has to stretch 
if the dragging force F is further increased. The satellite 
d-cone only occurs when the stretching component of the 
force becomes comparable with the bending component. 
Matching these two forces yields n = 3, in fair agreement 
with the in-silico value. 

V. CONCLUDING REMARKS 

We have shown that under generic conditions, the me- 
chanical response of an elastic membrane exhibits a sin- 



gular evolution. This response is a consequence of bend- 
ing minimization, the frictional collision between folds, 
and the birth of conical singularities. 

In particular, we identified the satellite d-cone as a ba- 
sic deformation mode of a moderately folded membrane 
when further constrained, which can be effectively de- 
scribed as a conical membrane with a larger thickness 
determined by the number of the folded layers. This 
is further supported by our experiments starting with a 
single-layer conical membrane. Our preliminary results 
and theoretical analysis show that the location of the tip 
of the satellite d-cone depends on the opening angle of 
the cone. 

Based on our experimental observations and numeri- 
cal simulations, we suggest a new route to model the me- 
chanical response of crumpled surfaces, by considering an 
effective system that preserves the overall symmetries of 
the deformed one and that has a renormalized bending 
stiffness. It seems feasible to use this approach to ex- 
pedite numerical simulations of crumpling and to make 
easier theoretical analysis of crumpled surfaces. 

We also showed that for any real and in-silico mem- 
branes the stretching mode dominates until the stretch- 
ing energy is enough to produce a global change, leading 
to the birth of the d-cone. Before the formation of the d- 
cone we showed that there is a universal scaling function 
for the dragging force that shows an anomalous exponent 
in the vertical displacement. 

Finally we showed that friction does play an impor- 
tant role in modifying the mechanical response of the 
membrane. Its effects emerge soon after the first colli- 
sion event between the two branches of the d-cone. We 
showed that friction will generically enhance the above 
mentioned renormalization process. It is expected that, 
as friction is increased, a d-cone appears sooner. This 
was confirmed by our experiments on membranes coated 
with materials having different frictional properties. 

We expect our findings to be applicable to enhance 
energy dissipation in structures under extreme deforma- 
tions by suitable modification of its frictional properties. 
For example, in car collisions more energy could be dissi- 
pated by designing the inner surface of metallic sheets to 
be structured in such a way that as soon as metal plates 
touch each other they get stuck. It would be interesting 
to see if this type of ideas can help us design structures 
with enhanced response to impacts. 
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